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AMENDMENTS TO THE CLAIMS : 

This listing of claims will replace all prior versions, 
and listings, of claims in the application: 

LISTING OF CLAIMS : 

1. (currently amended) A computation method for 
simulating behavior of a flow interacting with an object, 
comprising the steps of: 

the method providing 

developing a simulated numerical representation in N 

dimensions, N > 3 , that is- composed of a plurality of 

approximated values at a multitude of points in at least a part 
of space where the a flow interacts with the an object, the 
approximated values being of a physical parameter u of the flow 
to which is associated a velocity field a which determines a 
preferential direction [[,]] by mcano of executing a numerical 

scheme wherein at least one spatial p th derivative D p , p>l, of 

the parameter u is approximated at the points of the part of 
space, through the method comprising the steps of-=- 

i) constructing a discrete N -dimensional grid 

constructed by N families of coordinate lines and using the 
constructed grid using for the part of space^ a — diocroto -fcf— 
dimensional grid constructed by - N -familics of coordinate lines; 
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ii) computing, in at least one point P of the grid, 
called the point of computation, an approximated value Dp of D p 
with an error e n , by using values u s of the parameter in a 
collection of grid points, called the stencil S, and 
computational functions, evaluated with the values u s , the 
computational functions depending on the numerical framework in 
which Dp is expressed, and 

iii) choosing the computational functions for the 
approximated value Dp in such a way that the approximated value 
Dp is optimized for the preferential direction, and 

wherein the stencil S contains at least one point 
situated outside all the coordinate lines passing through the 
point of computation P , and the stencil S contains at least a 
first point and a second point, the first point being defined by 
N first coordinate lines of the N families of lines, the second 
point being defined by N second coordinate lines of the N 
families of lines, and for at least one family Nf of the 
coordinate lines, the first coordinate line belonging to the 
family Nf is different from and not adjacent to the second 
coordinate line belonging to the same family Nf ; a*** 
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outputting the simulated numerical representation that 
simulates, for the part of space, the behavior of the flow 
interacting with the object ; and 

determining the behavior of the flow interacting with 
the object by using the output simulated numerical 
representation . 

2. (original) The method of claim 1, wherein for 
varying preferential directions, the approximation depends 
continuously on the values u s . 

3. (original) The method of claim 1, wherein the 
computing step of Dp comprises the steps of: 

providing a local basis B(?i,?2/?3/— ) of curvilinear 
coordinates which has the unit vector along the preferential 
direction, and 

choosing the computational functions so that a 
contribution to the error e n of at least one pure or one mixed 
derivative as expressed in the local basis B is minimized, while 
using as a formulation of the values u s of the parameter at each 

of the points of the stencil S, a truncated Taylor series 
expansion with respect to the point of computation P with an 
error, called the truncation error e s . 
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4. (original) The method of claim 3, wherein the 
computational functions are individual coefficients C s and the 

approximated value Dp is a linear combination of values u s of 

the parameter. 

5. (original) The method of claim 3, wherein the 
computing step of Dp comprises the steps of: 

using an integral formulation for the computation of 

the derivative D p with computational functions which are fluxes 

through a control volume, 

computing in at least one volume an approximated value 

Dp of D p with an error s n , where the approximated value is a 

function of the flux formulation employed. 

6. (original) The method of claim 3, wherein the 
computing step of Dp comprises the steps of: 

subdividing the part of space into elements, containing 
nodes at which the approximation of the physical parameter u is 
stored, 

defining basis computational functions q> , called 
interpolation functions, on the elements, where the interpolation 
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functions are used to approximate the physical value u on the 
element, 

computing the integral of the derivative Dp on the 
element, using a test computational function y/ , called weighting 
function, 

computing in at least one element an approximated value 

D^'^ of Dp with an error e n , where the approximated value is a 

function of the interpolation functions <p and the weighting 
function (// . 

7. (original) The method of claim 3, wherein the 

computing step of D p comprises the steps of: 

subdividing the part of space into elements, containing 
nodes at which the approximation of the physical parameter u is 
stored, 

defining basis computational functions <p , called 
interpolation functions, on the elements, where the interpolation 
functions are used to approximate the physical value u on the 
element, 

computing the integral l e \ of the derivative Dp on 

volumes, and distributing parts cc{l e i to nodes i, where OCj 
represent computational functions called distribution 
coefficients . 
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8. (currently amended) The method of claim 3, wherein 

the computing step of Dp comprises the steps of: 

subdividing the part of space into elements, containing 
nodes at which the approximation of the physical parameter u is 
stored, 

defining basis computational functions <p , called 
interpolation functions, on the elements, where the interpolation 
functions are used to approximate the physical value u on the 
element, 

computing fluxes / at the surfaces of volumes, and 
distributing parts a { f to nodes i, where a { represent 
computational functions called distribution coefficients^ 
t [, ] ] 

9. (original) The method of claim 1, wherein the 
computing step of comprises the steps of: 

using a representation of the numerical solutions in 
Fourier components, and 

choosing in the approximated value D p the 

computational functions in such a way that the Fourier components 
are optimized for certain directions which are related to the 



13 



Docket No. 0594-1001-1 
Appln. No. 10/779,685 

velocity a, while using the values u s of the parameter at each of 
the points of the stencil S in the Fourier representation. 

10. (original) The method of claim 9, wherein the 
computational functions are individual coefficients C5 and the 

approximated value Dp is a linear combination of values u s of 
the parameter. 

11. (original) The method of claim 9, wherein the 

A 

computing step of Dp comprises the steps of: 

using an integral formulation for the computation of 

the derivative Dp with computational functions which are fluxes 

through a control volume, 

computing in at least one volume an approximated value 

D^ of Dp with an error e n , where the approximated value is a 

function of the flux formulation employed. 

12. (original) The method of claim 9, wherein the 
computing step of Dp comprises the steps of: 

subdividing the part of space into elements, containing 
nodes at which the approximation of the physical parameter u is 
stored, 
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defining basis computational functions <p , called 
interpolation functions, on the elements, where the interpolation 
functions are used to approximate the physical value u on the 
element, 

computing the integral of the derivative D p on the 
element, using a test computational function y/ , called weighting 
function, 

computing in at least one element an approximated value 
of D p with an error e n , where the approximated value is a 
function of the interpolation functions <p and the weighting 
function if/ . 

13. (original) The method of claim 9, wherein the 
computing step of D p comprises the steps of: 

subdividing the part of space into elements, containing 
nodes at which the approximation of the physical parameter u is 
stored, 

defining basis computational functions <p , called 
interpolation functions, on the elements, where the interpolation 
functions are used to approximate the physical value u on the 
element, 

computing the integral I d of the derivative D p on 
volumes, and distributing parts ajl e i to nodes i, where a { 
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represent computational functions called distribution 
coefficients . 

14. (currently amended) The method of claim 9, wherein 
the computing step of comprises the steps of: 

subdividing the part of space into elements, containing 
nodes at which the approximation of the physical parameter u is 
stored, 

defining basis computational functions q> , called 
interpolation functions, on the elements, where the interpolation 
functions are used to approximate the physical value u on the 
element, 

computing fluxes / at the surfaces of volumes, and 
distributing parts a { f to nodes i, where a { represent 
computational functions called distribution coefficients^ [[,]] 

15. (original) The method of claim 1, wherein the N- 
dimensional grid is expressed in a coordinate system which is 
chosen from the group consisting of: rectangular coordinates, 
spherical coordinates, cylindrical coordinates, parabolic 
cylindrical coordinates, paraboloidal coordinates, elliptic 
cylindrical coordinates, prolate spheroidal coordinates, oblate 
spheroidal coordinates, bipolar coordinates, toroidal 
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coordinates , 



conical 



coordinates, 



conf ocal 



ellipsoidal 



coordinates and conf ocal paraboloidal coordinates. 



16. (original) The method of claim 1, wherein the N- 



dimensional grid is chosen from the group consisting of: a grid 
with non-uniform mesh spacing, a grid which is moving, a grid 
which is deforming, a grid which is rotating, and a grid which is 
staggered, and any combination thereof. 

17. (original) The method of claim 1, wherein the 
spatial p th derivative D p is a pure derivative . 




18. (original) The method of claim 1, wherein p = \ - 



19. (original) The method of claim 1, wherein p = 2. 



20. (original) The method of claim 1, wherein the 



is a mixed derivative 



- with 



spatial p derivative D 



p 1 +p 2 +- = P • 
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21. (original) The method of claim 3, computing an 
approximated value D^of D 1 , where the approximated value is 
denoted by 



(*y)i,;,fc,... 
(u z )i,j,k,... 

V : 

wherein 



(27) 



( u x )i,j,k,- t*-m u i-m,j,k,. . - + a -m+l u i-m+l,i,k,. . . + • ■ • 

+ a n-l u i+n-\,jX--- +a n u i+n,j,k,... + 

( u y)i,j,k,... = ~{*-m( u i-m,j-m,k,... ~ u i-m f j f k f ...) 

+ a -m+l( u i-m+l,j-m+l,k,... ~ u i-m+l, j,k, + - - 
+ a n -l( u i+n-l,j+n-l,k,.-- - u i+n-l,j,k,...) 
+ a n( u i+n,j+n,k,... - u i+n,j,k,...)- hT y}' 

( u z)i,j,k,... = —^-m( u i-m,j-m,k-m,„. - u i-m,j-m,k,...) 

+ a -m+l( u i-m+l,j-m+l,k-m+l,... ~ u i-m+l,j-m+l,k t :.) + 
+ a n-l( u i+n-l,j+n-l,k+n-l,... - u i+n-l,j+n-l,k,...) 
+ a n ( u i+n,j+n,k+n,... ~ u i+n,j+n,k,...) + T z }, 

1 (28) 

wherein a_ m * 0 , fl n *0, m and n are given integers, 
m + n>0 , and m + n>l if m*n=0, the terms T x ,T y ,T z ,... represent 
the degrees of freedom which are used in the optimization of the 
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approximated value , and where indices (i,j f k,...) define the 

point of computation P on the N -dimensional grid, and 
Ax,Ay,Az,... denote the mesh spacings of the N -dimensional grid in 
each coordinate direction. 

22. (original) The method of claim 21, computing an 

approximated value of D x = — by a discretization of order 

de x 

d M+1 u 

M , in which the terms — u u with M 2 +M 3 + ... = M + l are 

de^de™*-- 

eliminated in the approximated value . 

23. (original) The method of claim 21, computing an 
approximated value of A = g^~ b V a discretization of order 

a M+1 w 

M , in which the terms — ^ M M with M 1 + M 2 + M 3 + . . . = M + 1 

and M 1 <M + 1 are eliminated in the approximated value D 1 in the 
case that is along the x-axis or along diagonals. 

24. (original) The method of claim 1, computing an 
approximated value Df of D 1 , where the approximated value is 
denoted by 
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("y )*,/,*,.. . 
(u z )i,j,k,... 



wherein 



(29) 



K )i,j,k,~ - = ^ Wm u i-m,j,k t . • • + *-m+l u i-m+l,j,k,. • . + • ■ ■ 
+ On-lUi+n-ltjX. • • + a n u i+n,j,k,-. }' 



(Uy)i,j,k,... =-^{t-rn(Ui-m,j-m,k,--- - u i-m,j,k,...) 

+ • a -m+l( u i-m+l,j-m+l,k,~. ~ u i-m+l,j ,k,...) + *• 
+ 0 n -l( M i+n-l,7+n-l,*,... - M i+n-l,j,k,...) 
+ M M i+n,;+n,k,... - M i+n,/,*,...)}' 



( u z)i,j,k,... ~ —7P-m( u i-m,j-tn,k-tn,... u i-m,j-m,k,...) 

+ a -m+l( u i-m+l,j-m+l,k-m+l,... - u i-m+l,j-tn+l,k,...) + ••• 
+ a n _ 1 (M i+n _ 1/7+n _i / ^ +n _i / ... -ttj+n-l,/+w-l,k,...) 
+ a n( u i + n t j+n,k+n,... ~ u i+n,j+n,k,...)\' 

: (30) 
wherein a_ m * 0 , 0 n *O, m and n are given integers, 
m + n>0, and ra + n>l if m*n=0, and where indices (i,j,k,...) 
define the point of computation P on the N -dimensional grid, 
and Ax,At/,Az,... denote the mesh spacings of the N -dimensional 
grid in each coordinate direction. 
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25. (original) The method of claim 21, wherein 
("* )«,y,k,. . . = - - - "i- W,*,. • • ) +T * } ' 

(u y ) irj ,k,... = ^{f ("i + l,y + l,k,.. + "i-l,y,k,... -"i-i,;-i,fc,..) +T y}' 

1 fl 

(u z )i,j,k,... ^—^ Ui+1 'j +1 ' k+1 ' - ~ Ui+1 ' j+1 ' k '- - 

+ Ui-\,j-l,k,... - "i-l,;-l,fc-l,.. .) + T z }' 

; (3D 



26. (original) The method of claim 21, wherein 
(%)i,y,k,...=^{^("»-2^^ 
1 f 1 

("y)t,;,fc,... = ^{^ (Ul - 2 '/'- 2 ^ - - M «'-2,;^ - 
- 8u f _ 1/; _i /fc/ . . . + 8ui_ h j ik> _ . . 
+ 8uj+i,y+i,k,. . . - 8ui+l,j,k,. . . 
~ "i+2,/+2,k,. . . + u i+2,j,k,. ■ - ) + T y }' 

1 f 1 

(" Z )i,;,k,... = — j^("i-2,y-2,fc-2,... -"•-2 / ;-2,k,... 

- 8«i-i,y-i,k-i,. . . + 8 "f-i,y-i,fc,. . . 

+ 8u,+i,y+i,k+l,. . . - 8"i+i,y+l,k,. . . 

- "i+2,y+2,k+2,. . . + u i+2,j+2,k,. .) + T z)> 

! (32) 



27. (original) The method of claim 21, wherein 
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1 f 3 1 



("y )*,;,*,... =^{- 2 ( u i-lJ-l,k,...-Ui-l,j,k,.*.) 

+ ^( u i-2J-2,k,...-Ui_ 2 JX-.> T y}> 



( u z )i,j,k,. . . = ^ {" 2 ( u i-l,j-l,k-l,. . . ~ . . ) 

+ ^( w i-2,;-2,*-2,... -«i-2,;-2,*,.. .) + T z} ' 



(33) 



28. (currently amended) A oimulation method according 
to claim 1, in three dimensions, obtaining the approximation 

with order M of the derivative D v = with 

Vl + Pl + P3~P on a <? r id of given extent from the — output — — fche 
program gcncratc - discrctizations which io given in appendices 1 - 5 
from the output of a computer program. 

29. (currently amended) A oimulation method according 
to claim 28, wherein the input parameter optimize has the value 1 

or 2 in the approximation by a discretization of order M 
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o u 

terms are eliminated, which are the terms — — 0 with 

de™ 2 de™* 

d M+l u 

M 2 + M 3 + ... = M + 1 , or the terms terms TmTTmTTW, — with 



M 1 +M 2 +M 3 +... = M + 1 and M 1 <M + 1 in the case that i?i is along 
the x -axis or along diagonals. 

30. (currently amended) A simulation method according 
to claim 28, wherein the input parameter order hao the value 1 

the approximation of is a first order accurate 

discretization . 



31. (original) The method of claim 3, computing in 
three dimensions an approximation of the second derivative 

D 2 = , where the approximation D 2 is expressed in the terms 

del 

u xx ,u yy ,u zz ,u xy ,u yz and u zx which are given by 

("xx)i,j,k =—72 ("i+W,* - 2u i.j,k +u i-l,j,k +T xx) > 
(Ax) 

\ 

( u yy)i,j,k = -—T(ui,j+l,k ~ 2u i,j,k + +T yy) > 

) 

i 

(Uzz)i,j,k =-—2(Ui,j,k+l- 2u i,j,k +u i,j,k-l +T zz) - 
(Az) 
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1 fl 

(u>xy)i,j,k = ^^{2^ Mf+1 ' ; ' +1 '* ~ Ui >i +1 > h - U i+l,j*+ 2u i,j>k 
- u i-l,j,k - u i,j-l,k+ u i-l,j-l,k) + T xy}> 

1 Jl 

i u yz)i,j,k "^^j^+V+l^+l ~ u i+l,j,k+l - u i+l,j+l,k + u i+l,j,k 

+ M i,;+1,*+1 ~ u i f jMl ~ u i,j+l,k+ 2u i,j,k 
~ u i,j-l,k ~ u i,j,k-l+ u i,j-l,k-l 

+ -Ui_ X j fk _ x +Ui_ lf j tk ) + T yz }, 



1 (1 

( u zx)i,j,k = '£^\j( u i+l.j+l.k+l ~ u i,j+l,k+l - u i+l,j+l,k+ u i,j+l,k 
+ -u i+1/ j /k -u if j M1 +2u it j fk 

- U i-l,j,k ~ Ui,j,k-l + M i-l,/,Jt-l 

+ Mf_i /; _i /fc „i - Ui_ lf j_ lfk - u it ^ x k _ x + WiJ-l,* ) + T zx } , 

(34) 



wherein the terms T xx ,T xy ,T xz ,Tyy,T yz and T 2Z 
represent the degrees of freedom which are used in the 
optimization of the approximated value , and where indices 

(i,j,k) define the point of computation P on the three-dimensional 
grid, and A*,Ay,Az denote the mesh spacings of the three- 
dimensional grid in each coordinate direction. 



32. (original) The method according to claim 1, wherein 
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Dp = J^ m L n Dp n where L n are constants and each Dp n is 
a function of values u s of the parameter in a collection of grid 
points, called stencil S n , with individual computation 

functions, which depend on the numerical framework in which Dp n 

A 

is expressed, and wherein in the approximation Dp n , the 
computation functions are chosen in such a way that the 
approximated value Dp U is optimized for the preferential 
direction . 

33. (original) The method of claim 1, wherein 

A , A 

Dp = / in L n Dp n where L n are limiting functions of the 

values u s of the stencil S, and at least one Dp n is a function 
of values u s of the parameter in a collection of grid points, 
called stencil S n , with individual computation functions, which 

depend on the numerical framework in which Dp n is expressed, and 

wherein in the approximation Dp n , the computation functions are 

chosen in such a way that said approximated value Dp n is 
optimized for the preferential direction. 

34. (original) The method of claim 1, wherein the 
stencil S is chosen from the group consisting of: upwind 
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discretization stencils, centered discretization stencils, and 
discretization stencils which are a combination of at least one 
upwind discretization stencil and at least one centered 
discretization stencil . 

35. (original) The method of claim 1, wherein the 
numerical discretization is a non-linear discretization. 

36. (original) The method of claim 1, wherein the 
numerical scheme is chosen from the group consisting of: the Lax- 
Wendroff scheme, the Lax-Friedrich scheme, the MacCormack scheme, 
the leap-frog scheme, the Crank-Nicholson scheme, the Stone-Brian 
scheme, the box scheme, Henn 1 s scheme, the QUICK scheme, the 
/rscheme, the Flux Corrected Transport (FTC) scheme, the family 
of ENO schemes, schemes in the class of the Piecewise Parabolic 
Method (PPM) , multi-level schemes, and schemes obtained with the 
fractional step method and variants thereof. 

37. (original) The method of claim 1, wherein the 
numerical scheme includes the discretization of a plurality of 
equations . 

38. (original) The method of claim 1 for the numerical 
simulation of physical phenomena which are modeled by the Navier- 
Stokes equations with equation (s) of state. 
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39. (original) The method of claim 1 for the numerical 
simulation of physical phenomena which are modeled by the Euler 
equations with equation (s) of state. 

40. (original) The method of claim 1 for the numerical 
simulation of physical phenomena which are modeled by the 
magneto-hydrodynamic equations with equation (s) of state. 

41. (original) The method of claim 1 in combination 
with at least one model for the physical phenomena chosen from 
the group consisting of: equations to model turbulence, equations 
to model chemical reactions, equations to model electromagnetic 
phenomena, equations to model multiphase flow and equations to 
model multi-physics phenomena. 

42. (original) The method of claim 1 for the numerical 
simulation of physical phenomena in combination with at least one 
acceleration technique chosen from the group consisting of: local 
time-stepping, multi-grid, GMRES and preconditioning. 

43. (original) The method of claim 1 for the simulation 
of a physical phenomenon which includes a material flow. 
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44. (original) The method of claim 1 for the simulation 
of a physical phenomenon which includes a material object 
interacting with material flow. 

45. (currently amended) The method of claim 1 , wherein 
the object is for the oimulation of a phyoical phenomenon which 
includco a vehicle. 

46. (currently amended) The method of claim 1 , wherein 
the object — fe^e — oimulation — — a — phyoical — phenomenon — which 
includes a rotating blade. 

47. (currently amended) The method of claim 1 , wherein 
the object is — fehe — oimulation — e# at least a part of the 
atmosphere of . the earth. 

48. (currently amended) The method of claim 1 , wherein 
the flow includes oil for the numerical oimulation of . 

49. (currently amended) The method of claim 1 , wherein 
the flow includes #en? — fehe — numerical — oimulation — — phyoical 
phenomena which include combustion. 

50. (currently amended) A data processing system 
programmed to implement a oimulation method according to claim 1. 
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51. (original) A computer program that can be loaded in 
a data processing system so as to implement a method according to 
claim 1 . 

52. (currently amended) A digital storage computer- 
readable medium containing a stored computer program that is 
configured to be loadable in a data processing system so as to 
control the data processing system to implement a method 
according to claim 1 . 
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